Radio Cherenkov signals from the Moon: neutrinos and cosmic rays 
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Neutrino production of radio Cherenkov signals in the Moon is the object of radio telescope 
observations. Depending on the energy range and detection parameters, the dominant contribution 
to the neutrino signal may come from interactions of the neutrino on the Moon facing the telescope, 
rather than neutrinos that have traversed a portion of the Moon. Using the approximate analytic 
expression of the effective lunar aperture from a recent paper by Gayley, Mutel and Jaeger, we 
evaluate the background from cosmic ray interactions in the lunar regolith. We also consider the 
modifications to the effective lunar aperture from generic non-standard model neutrino interactions. 
A background to neutrino signals are radio Cherenkov signals from cosmic ray interactions. For 
cosmogenic neutrino fluxes, neutrino signals will be difficult to observe because of low neutrino flux 
at the high energy end and large cosmic ray background in the lower energy range considered here. 
We show that lunar radio detection of neutrino interactions is best suited to constrain or measure 
neutrinos from astrophysical sources and probe non-standard neutrino-nucleon interactions such as 
microscopic black hole production. 



I. INTRODUCTION 

Ultrahigh energy neutrinos may originate from astro- 
physical sources, from exotic sources such as ultramassive 
particles which decay and from cosmic ray interactions 
with the background radiation [T]. Because cosmic rays 
have been observed up to energies of 10 19 eV, high energy 
neutrino flux from cosmic ray interactions with photons 
producing charged pions [2] which decay into neutrinos 
[3] is a "garanteed" neutrino flux. Deatiled flux predic- 
tions of these final state neutrinos, called GZK-neutrinos 
or cosmogenic neutrinos, still have theoretical uncertain- 
ties associated with the composition and injection spec- 
trum of the highest energy cosmic rays, and the photon 
spectum now and at earlier epochs, as discussed in, e.g., 
Refs. 4-9 . Observations of neutrino signals will be an 
important piece of the high energy astrophysics and cos- 
mological picture. 

In contrast to photons, neutrinos have weak interaction 
cross sections [I014T4] so neutrino fluxes are not attenu- 
ated over cosmic distances. However, one needs detectors 
sensitive to many targets for neutrino interactions since 
the interaction probability is low. There are already a 
number of effort to observe these highest energy neutri- 
nos, including signals from the air showers they would 
produce [T5], and from the particles and radiation they 
produce when they interact in matter [T6Tf2"5] . All obser- 
vational efforts require large volumes. 

Among these observational efforts are neutrino induced 
events on the Moon, where the Moon is the target, 
and the signal is the radio Cherenkov emission [T9Tl25| . 
Neutrinos, when they interact with nucleons and nu- 
clei, generate hadronic showers with shower energies of 
_E s h r ~ 0.2_E„. An electron charge excess is produced, 



and this group of electrons moves faster than the speed 
of light in the lunar regolith, hence the Cherenkov signal 
[26]. 

Evaluation of the signal depends on the cosmogenic 
neutrino flux, the ultrahigh energy neutrino cross section, 
the radio signal production, attenuation and refraction at 
the surface and detection parameters. An approximate 
expression for the effective aperture for neutrino induced 
radio Cherenkov signals from the Moon has been devel- 
oped by Gayley, Mutel and Jaeger(GMJ) in Ref. [27] . 
Ultrahigh energy neutrinos are incident isotropically on 
the Moon to a good approximation. In analogy with 
terrestrial observations, e.g., at the IceCube neutrino ob- 
servatory [28] . neutrinos that are incident on the "back- 
side" of the Moon and traverse a portion of the Moon 
before interacting are called "upward" neutrinos. Neu- 
trinos which interact on the surface of the Moon facing 
the Earth are denoted "downward" neutrinos. 

Gayley et al. find that, depending on the energy and 
cross section, the event rate from neutrinos incident on 
the surface we see (downward neutrinos) sometimes dom- 
inate over the event rate from neutrinos incident on the 
surface of the Moon not visible from Earth (upward neu- 
trinos) 27 . This may seem counter-intuitive, however, 
between the various angles of incidence, the Cherenkov 
angle and the angular spread of the Cherenkov cone, an- 
gles of refraction and angle characterizing the lunar sur- 
face roughness, the downward neutrinos can contribute 
appreciably to the signal. 

Cosmic ray interactions in the lunar regolith produce 
hadronic showers as well. These cosmic ray induced 
hadronic showers also induce a charge excess and a re- 
lated radio Cherenkov signal which can be used to mea- 
sure or constrain the ultrahigh energy cosmic ray flux 
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FIG. 1: The flux of cosmic rays (solid) with the parametriza- 
tion from Ref. [3U] and three models for the flux of cosmo- 
genic neutrinos (dashed) from Ref. [9] with the assumption 
of proton primaries. 



[29] , The cosmic ray flux is quite large at low energy, 
but it falls off rapidly with increasing energy. The flux 
of cosmic rays (solid line) and three models for the flux 
of cosmogenic neutrinos from Ref. [9] (dashed lines) are 
shown in Fig. [TJ The three models for the cosmogenic 
neutrinos are a subset of predictions for the cosmogenic 
neutrino flux. These three come from the initial assump- 
tion of cosmic rays which are primarily protons, in which 
the sources are assumed to have maximum acceleration 
energies of E p ^ max = 10 11 GeV, 10 115 GeV and 10 12 
GeV. Details of the model, including injection spectra, 
assumptions about the star formation rate and the cor- 
related UHE cosmic ray sources and other inputs to the 
evaluation are described in Ref. [S]. The cosmic ray flux 
shown uses the parametrization of the Auger Collabora- 
tion in Ref. [15]. 

Given the cosmic ray flux |30) and its production of ra- 
dio Cherenkov signals similar to the corresponding neu- 
trino Cherenkov signals [2§], we evaluate the cosmic ray 
induced signals. The considerations for cosmic rays, as 
compared to neutrinos, are somewhat different. Cosmic 
ray fluxes are attenuated in the Earth's atmosphere, how- 
ever, there is no lunar atmosphere. Within the lunar 
regolith, the cosmic ray interaction length is small com- 
pared to any other characteristic distance scale. There 
will not be any "upward" cosmic ray induced signals, 
but there is the potential for "downward" cosmic rays to 
produce radio Cherenkov signals. Cosmic ray flux atten- 
uation in the lunar regolith is, of course, a large effect, 
resulting in cosmic ray signals being produced very close 
to the lunar surface. In Ref. [33], ter Veen et al. have 
determined that the radio Cherenkov signal far from the 
lunar surface is essentially the same as one originating 
deeper in the regolith (apart from attenuation). 

In this paper, we look at the relative importance of 



cosmic ray and neutrino induced radio Cherenkov signals 
from the Moon assuming standard model cross sections 
and the neutrino and cosmic ray fluxes in Fig. 1. We use 
the approximate expression from Ref. |27j for the event 
rates, and we modify this analytic result to account for 
strong attenuation where applicable. We also consider 
non-standard model neutrino-nucleon cross sections and 
generic neutrino fluxes from astrophysical sources. 

There have been a number of discussions of the cross 
section dependence of various neutrino induced signals 
|31H33j . In the next section, we review the neutrino cross 
section and the resulting dependence of the event rate on 
the cross section. In Section III, we evaluate the event 
rates using the standard model neutrino nucleon cross 
section and for some rescaled cross sections, all using the 
approximate analytic GMJ expression of Ref. [27 • The 
dependence of the event rate on detection characteris- 
tics, namely the radio frequency and minimum detectable 
electric field, is shown. 

In Section IV, we show how the approximate analytic 
expression is modified to account for cosmic ray flux at- 
tenuation on short distance scales in the lunar regolith. 
We find that the event rate is independent of cosmic ray 
cross section, as long as the cosmic ray interaction length 
is small compared to the photon attenuation length. We 
compare the rates of cosmic ray and neutrino induced 
radio signals using the fluxes shown in Fig. 1 and the 
standard model neutrino nucleon cross section. We show 
in Section V how the inclusion of mini-black hole produc- 
tion, as an example of a neutrino-nucleon cross section 
enhancement, and how alternative neutrino spectra affect 
the predicted event rates. 

Our conclusions appear in Section VI. We find that 
while lowering the minimum electric field detectable by 
a radio telescope array would help increase the number 
of neutrino events, since it effectively lowers the neutrino 
energy threshold for detection, it also increases the num- 
ber of cosmic ray events. Cosmogenic neutrino fluxes on 
the scale of those presented in Ref. [S] will be difficult to 
observe on the one hand because of low fluxes (at the high 
energy end) or because of the cosmic ray background (in 
the lower energy range considered here) . Lunar radio de- 
tection of neutrino interactions is best suited to constrain 
or measure neutrino sources other than the cosmogenic 
sources and non-standard neutrino nucleon cross section 
enhancements. 



II. NEUTRINO CROSS SECTIONS AND 
EFFECTIVE SOLID ANGLE 

The dependence of neutrino induced events on the neu- 
trino cross section has been the subject of much dis- 
cussion, e.g., in Refs. [3TH33] . One feature is that the 
probability of interaction is proportional to the neutrino 
nucleon cross section, however, the neutrino flux attenu- 
ation is also affected by the cross section. 

Experimentally, the neutrino nucleon cross section has 
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been directly measured for E v < 450 GeV [Mj. A re- 
lated measurement, the charged current interaction cross 
section for electrons in ep scattering at HERA translates 
to a neutrino cross section with E v = 27 TeV incident 
on a proton at rest [35] . The moment transfer relevant 
to UHE neutrino scattering is characterized by the W- 
boson mass M^y. At Q 2 ~ M™-, the structure functions 
have been measured in the Bjorken x regime of x larger 
than a few times 10~ 3 . At the Large Hadron Collider, 
one expects measurements of the structure functions for 
x > 10~ 5 for similar values of Q 2 . The reach in x at the 
Large Hadron Collider (LHC) extends to an equivalent 
neutrino energy E v ~ 10 8 GeV, using the approximate 
correspondence that 2M N E V ~ ~ (Q 2 ) Wc 

consider here neutrino energies E v > 10 8 GeV, above 
kinematic regions probed even by the LHC experiments. 

The neutrino cross section in the standard model, 
with a power law extrapolation of the structure 
functions at low x, has an approximate power 
law energy dependence [TP! [TT] cr*^ ~ 1.57 x 
10~ 35 cm 2 (_E I/ /GeV) - 33 . The neutrino interaction length 
is defined to be 
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FIG. 2: The neutrino interaction length (in centimeters water 
equivalent distance) as a function of energy, with indications 
of the diameter of the Earth and Moon using the cross sec- 
tions from Ref. |12| . 
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This has the dimensions of g/cm 2 , or equivalently cen- 
timeters of water equivalent distance (cmwe), which can 
be converted to a distance by dividing by density. Rather 
than use a power law form, we use the charged current 
cross sections of Ref. [IS], with the addition of the top 
quark contribution (on the order of a 30% correction at 
the highest energy [l2j US]). The charged current cross 
section, for high energies, is multiplied by an approxi- 
mately energy independent constant of 1.43 to obtain the 
total neutrino-nucleon cross section [14]. The standard 
model neutrino interaction length is shown in Fig. [2] For 
an incident neutrino energy of 10 3 GeV, the tau neutrino 
cross section is only ~ 5% less than the muon neutrino 
cross sections [36] . We will be considering higher ener- 
gies where we can use the same cross sections for all three 
neutrino flavors. 

The importance of neutrino attenuation in the Moon 
or Earth can be seen with a comparison of the interac- 
tion length L v as a function of neutrino energy for the 
standard model cross section with the Earth and Moon 
diameters. The Earth diameter equals the neutrino in- 
teraction length for E v larger than a few tens of TeV. 
The Moon's diameter (in cmwe) equals the neutrino in- 
teraction length for an incident neutrino energy of a few 
times 10 6 GeV. Our focus will be on cosmogenic neutri- 
nos at energies much higher than 10 6 GeV, so attenu- 
ation effects will be important in the evaluation of in- 
teraction rates, as we discuss in detail below. We note 
that the Earth's horizontal atmospheric depth is small 
on the scale of the neutrino interaction length, even at 
E u = 10 12 GeV. 

The cross section dependence of event rates is specific 



to the signal and the possibility of neutrino regeneration 
through neutrino neutral current interactions and for tau 
neutrinos, through neutrino production and decay [311 
1321 137[ . The neutral current neutrino regeneration can 
be approximately included via an effective cross section 



= a vN 



(2) 



where we have taken a to scale with energy as cr°^. In 
principle, k depends on neutrino energy, but in practice 
at high energy, the differential cross section for neutrino 
neutral current interactions has an approximate scaling 
with neutrino energy the same way as the total cross 
section. We neglect the v T — > r — » v T regeneration, which 
is typically not very important for steeply falling fluxes 

m\- 

We can estimate n in the standard model with some 
approximations. This is also outlined in the Appendix 
of Ref. [27] . The neutrino flux as a function of column 
depth X is 



dX 



$ v (E,X(0)) 

dE'$ v (E',X{9)) 



(3) 

da NC (E',E) 
dE 



In the standard model, the neutral current cross section 
is approximately r^vc = 0-3 of the total cross section. 
Using 



da NC (E',E) 
dE 



r NC 4t(E')S(E-(l-{y))E') (4) 



where y = (E' — E)/E' is the neutrino inelasticity, relat- 
ing the change in neutrino energy when it interacts, nor- 
malized to the initial neutrino energy. At high energies, 
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(y) ~ 0.2. For a power law spectrum $> V (E, X{9)) ~ _E 7 
and for a neutrino cross section which scales with energy 



as a 



E , Eq. (Wb can be written as 



d$ v (E,X(6)) 
dX 



■(l-r NC (l-(y)) r '- 1 - s ) 
*v(E,X{fi)) 



$> V (E,X{9)) 



L v 



(5) 



When (5 = 0.3, as is approximately the case for the stan- 
dard model, with tmc — 0.3, n ~ 0.74 for a spectral 
index 7 = 2. The value of k is not very sensitive to the 
spectral index. It increases to k ~ 0.84 when 7 = 4. We 
use k = 0.84 in our evaluation below, and we define the 
attenuation distance A to include the regeneration effect 
by 



A = 1/(jN a = L v /k 



(6) 



We review here the scaling of upward event rates in- 
cluding attenuation as a function of the cross section. 
Schematically, event rates T for a detector of cross sec- 
tional area A are given by 

r= I dE v dQ. v dA-n{9 v )drV{E v ,9 v ,r) 

x § V {E V ,X(9 V )) (7) 

where <& V (E V ,X) is the neutrino flux in units of 
neutrinos/ (cm 2 s sr GeV), and drV(E„,9 v ) is the proba- 
bility the neutrino of energy E v produces a signal in the 
interval dr. The angle 9 V is the incident angle of the neu- 
trino flux with respect to vector normal to the cross sec- 
tional area of the detector. The probability to produce 
a signal depends linearly on the neutrino cross section 
for the specific signal (cr s ), so considering short distances 
L = J dr, we can write drV(E^,9^) — dra s NAP- The 
effective volume of the detector is V, where dV — drdA. 
We consider detection near the surface, where for upward 
neutrino fluxes, the depth of the detector is negligible. 

We start with a configuration of a detector near the 
surface which is approximately isotropic: where A ■ n 
is independent of the neutrino direction and where the 
pathlength of the neutrino in the detector of size V = L 3 
is approximately independent of incident neutrino direc- 
tion: 

V iso (E u ,9„,L) ~V iso (E v ,0,L) . 

The neutrino flux accounting for attenuation, for an up- 
ward neutrino traversing a sphere and emerging with an 
angle 9 V with respect to the normal to the surface is ap- 
proximately 



* v (E vt X{9 v )) ~ e- aJlcOBfl "'/ A * I ,(£„,0) 



(8) 



for a flux incident on the Earth or Moon with column 
depth of the diameter of 2R. For the Moon, 2R = duPu 
for the diameter of the Moon d^ j — 3, 480 km and average 



lunar density pu — 3.34 g/cm 3 . Here, 9 V is the zenith 
angle of the incident neutrino flux. We assume that the 
neutrino flux is isotropically incident on a spherical body 
(the Earth or Moon). 

By factorizing the depth dependent neutrino flux as in 
Eq. ^ , the integrand for an isotropic incident flux has a 
factor which includes the effective solid angle fl e s where 



O lso 
2tt 



7T/2 



d9 v sin 9 v e 



-2i?.cos e v /x 



2R\ 



(l - exp(-2i?/A) 



(9) 



The angular integral for a fixed energy (fixed A) de- 
pends only on 2R/X. This is shown in Fig. [3] where 
the dashed line shows the quantity X/2R. The dashed 
line is the scaling behavior of the integral for small inter- 
action lengths compared to the diameter of the Earth or 
Moon (X/2R < 1/2). 

For detectors which are not isotropic, similar result 
is found. Consider the a detector schematically shown 
in Fig. |4j where the area is A = A and the detector 
thickness is d, the quantity A ■ n — A cos 9 V appears in 
the integral, but now f drV(E v ,9 v ,r) — a s NApd/ cos#„, 
assuming no neutrino attenuation through the depth of 
the detector. The combined angular dependence makes 
the integrand the same as for an isotropic detector, so 
fijffj is the same for any isotropic upward neutrino flux 
incident on an underground detector. Eq. ([9| shows 
that roughly, for a given X/R, the zenith angles which 
contribute range between 9 = 90° and 9 = cos _1 (A/2i?) 
when 2R > A. 

Combined with the neutrino interaction probability 
and the effective area A = Aq, the cross section depen- 
dence of the effective solid angle, for small interaction 
lengths relative to the column depth 2R is 



drSy£P(E v ,0,r)A~Ao 



<J S dp 
~^2R 



(10) 



for isotropic incident fluxes. This agrees with the dis- 
cussion of, e.g., Ref. [31], where it is noted that for large 
neutrino cross sections, the effective solid angle is reduced 
in just the proportion to the increase in event rate in 
the detection region in their discussion of upward shower 
events. 

The slab configuration of Fig. [4] is relevant to neutrino 
induced radio Cherenkov signals. As we discuss below, 
radio signals are produced in the a thin layer of the lu- 
nar regolith on the surface facing the Earth. The depth 
of layer is characterized by the attenuation length of the 
radio signal, L 7 . Additional factors are required to ac- 
count for the radio signal production and refraction at 
the surface, but schematically, the upward signal will be 
independent of the neutrino cross section at high ener- 
gies, and it is proportional to L 7 . 
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FIG. 3: The effective solid angle divided by 2n for an under- 
ground detector with an isotropic effective ill Oil, OS £L function 
of the ratio of the diameter (in water equivalent distance) to 
the attenuation distance (solid line). The dashed line shows 
(X/2R). 




combines the neutrino interaction probability, effective 
area (A) and effective solid angle, 



r = J dE v A M {E v )$ v {E v ,X) . 



(11) 



Eq. (11) also accounts for the radio Cherenkov signal 



production, attenuation and refractions at the lunar sur- 
face. Am(E) has units of sr-km 2 . GMJ[27] have shown 
that it is possible to parametrize the effective aperture 
for a lunar radio Cerenkov signal with the analytic form 



A M {E) = A 



{nl - 1) L„ 



8n r L 
x (* ds + + * 

= Ads + Adr + A u 

= A P(E) . 



/o 3 A 



(12) 



In this equation for the Moon, the index of refrac- 
tion is n r = 1.73 and the maximum lunar aperture 
is Ao = Att(ttR 2 m ) for lunar radius Rm — 1,740 km. 
The photon interaction length (in g/cm 2 units) is L 7 = 
9 m x pj (iy/GHz) with a lunar regolith density of p = 1.7 
g/cm 3 . We discuss our neutrino interaction length in de- 
tail below, but this factor of L~ x carries the factor of a s 
in the probability for the neutrino to produce a signal. It 
is convenient to scale out a factor of A and look at the 
energy dependent function P(E). 

The parameter fo is the ratio of the thickness of the 
Cherenkov cone at the electric field threshold £ m ; n to the 
full thickness of the Cherenkov cone (2Aq), 



/o 
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0.6e 



(13) 



where already the requirement that the electric field of 
the signal at the Earth is larger than the electric field 
threshold of the detector has been enforced. The quantity 
So depends on the distance to the moon d = 3.84 x 10 5 
km, the energy of the neutrino induced shower E s h r and 
the radio frequency of the radiation v is 



FIG. 4: A detector of area A and depth d, with neutrinos 
incident in a directions characterized by n which makes an 
angle 9 V to the normal to the surface A. 



III. LUNAR RADIO CHERENKOV SIGNALS 
FROM NEUTRINOS 

To evaluate the neutrino event rate, we look at the ef- 
fective aperture of the Moon evaluated by Gayley, Mutel 
and Jaeger in Ref. [27|. The effective aperture Am{E u ) 



£ = 0.0845 



V 



mMHz 



2.32 GHZ 



EeV 
-l 



GHz 



(14) 



where -E s hr — 0.2E 1 ,, is the approximate hadronic shower 
energy for neutrino interactions with nucleons in the lu- 
nar regolith. The quantity Ao is the Cherenkov cone half 
width, 



A = 0.05 



GHz 



1 + 0.075 log 



E f 



hi 



10 



10 10 GcV 



(15) 
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The electric field threshold e m ; n of the detector de- 
pends on the collection area of the telescopes, the band 
width and other features of the specific telescope or ar- 
ray of telescopes [37]. One example is Project RESUN 
[2"T] using the radio Expanded Very Large Array where 
£min ~ 10~ 8 V/m/MHz at 1.4 GHz. 

There are three terms in eq. ( 12 ) representing the 



angular aperture fractions for downward neutrinos on 
a smooth surface (ds), downward neutrinos on a rough 
surface (dr) and upward neutrinos (u). They are respec- 
tively, 



/oA 
16 

3-372*0 = tan" \0.Uu - 22 ) 



16 



A 



16 L, 



3 2RmPm 3 2RmPmk 



(16) 



These terms were derived in the approximation that the 
neutrino interactions occur near the lunar surface, and 
that the spread of the Cherenkov cone is small, as is the 
incident neutrino direction relative to the horizontal to 
the lunar surface. This last approximation is valid as long 
as the neutrino interaction length is small compared to 
the lunar diameter. 

Qualitatively, the prefactor accounts for the incident 
angles of the neutrinos, the cross sectional area of the 
moon, the interaction probability in the outer layer of 
the lunar regolith 



P. 



eff 



^sinOJ 2 



the integral over the width of the Cherenkov cone giving 
sin 6> c /oAo and the refraction effect on the solid angle sub- 
tended by the emerging radiation from inside the Moon 
(n r ). We recall that sin 2 8 C — (n 2 — l)/n 2 r . 

For the upward effective aperture, the factor & u ac- 
counts for the reduced effective solid angle due to neu- 
trino attenuation in the Moon. In fact, 



16 fig 
3 2tt 



(17) 



If the neutrino nucleon cross section is smaller than the 
standard model cross section, or at lower energies, the 
scaling of fijfg with (X/2R) is modified, so we make the 
substitution of eq. ( Vty in the event rate to allow fol- 
lower cross sections (larger interaction lengths). 

For neutrinos that are incident downward, some radio 
signal will emerge, namely the portion of the solid angle 
equal to the thickness of the Cherenkov cone. Without 
surface roughness, this is the only contributor to the ra- 
dio signal, however, surface roughness permits the radio 
signal that would otherwise be lost to emerge. The ap- 
proximate analytic result of Gayley et al. has the ^dr 
contribution proportional to ctq = v2 tan- 1 (0.14i/ ' 22 ), 
the surface roughness parameter in terms of the radio 
frequency v in GHz. For v = 1.5 GHz, a = 12.3° = 0.21 



rad. 

In Fig. [5j we show the effective aperture for (a) 
v = 0.15 GHz and (b) v = 1.5 GHz, for e min = 10~ 8 
V/m/MHz. For the lower frequency, the downward neu- 
trinos in the smooth approximation dominate, while for 
the higher frequency, roughness on the lunar surface 
transmits radio signals that would otherwise be lost. At 
the lower frequency end, the angular spread Cherenkov 
cone is broader than at higher frequencies, allowing more 
of downward signal to head towards Earth, even with 
an approximately smooth lunar surface. For the higher 
frequencies, the downward signals that emerge from the 
smooth surface approximation (ds) are smaller than the 
downward signals accounting for surface roughness. This 
conies from /0A0 ~ 0.01 — 0.1 which is small compared 
to the solid angle characterizing the surface roughness, 
cto — 0.2, as discussed by Gayley et al. [27]. In Fig. |6j 
for e m ; n = 10~ u V/m/MHz, we show the effective aper- 
ture for the same two radio frequencies. The lower min- 
imum detectable electric field allows a probe of a lower 
neutrino energy, where the flux of neutrinos is predicted 
to be much larger. 

To show the dependence of the event rates on the neu- 
trino nucleon cross section, we take as an example the 
radio frequency of v = 1.5 GHz and 100 hrs of viewing 
time. In Figs. [7] and [8] we use a cross section scaling 
factor of S of the standard model (SM) neutrino nucleon 
cross section, 



q„SM 



(18) 



to evaluate the rates for two choices of £ m i n , e m i n = 10 8 
V/m/MHz and e min = lO" 11 V/m/MHz. For the higher 
electric field threshold, we have shown results only for 
the highest cosmogenic neutrino flux in Fig. [T] For the 
lower threshold, we show results for both the upper and 
lower cosmogenic flux predictions. 

Each of the separate contributions are shown: the dot- 
ted line show the "down smooth" contribution, the dot- 
dashed line shows the "down rough" contribution and 
the dashed line shows the "up" contribution. As dis- 
cussed, the "up" contribution becomes independent of 
neutrino nucleon cross section when the cross section is 
large enough, while the "down" contributions scale lin- 
early with the cross section. For low cross sections, atten- 
uation of the upward neutrino flux is less prominent, so, 
for example, with e min = 10~ n V/m/MHz and S ~ 0.01, 
the upward event rate scales with the neutrino nucleon 
cross section. 

The total numbers of events in one hundred hours are 
shown with the solid curves. The dots labeled CTW 
in each of the figures show a range of predictions using 
the uncertainty bands of Connolly, Thorne and Waters 
(CTW) discussed in Ref. [2]. The uncertainty is largest 
at the highest energies, probed with the lower electric 
field threshold. While the range of predictions spans a 
factor of about five, the overall predicted rate is quite 
low, even for the higher cosmogenic flux. The standard 
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FIG. 5: The effective aperture as a function of neutrino energy 
for e m i n = 10~ 8 V/m/MHz and for two frequency choices: (a) 
v = 0.15 GHz and (b) v = 1.5 GHz. The solid line shows 
the total of the downward rough (dr, dot-dashed), downward 
smooth (ds, dotted) and upward neutrino (up, dashed) con- 
tributions. 



FIG. 6: The same as Fig. 5 but for e min = 10 -11 V/m/MHz. 



field detection threshold at Earth must be on the order 
of e m i n = 10~ n — 10~ 10 V/m/MHz for even one event in 
100 hrs. 



model cross section (or the flux) must be enhanced by at 
least six orders of magnitude to get one predicted event 
in 100 hours. If it is the cross section that is enhanced 
to this degree, the neutrino interaction length becomes 
small compared to L 1 and the approximations used here 
for neutrinos do not apply. We discuss this possibility in 
the next section. 

If an electric field detection threshold can be as low 
as £min — 10~ n V/m/MHz, then on the order of one 
event is predicted for v = 1.5 GHz. The standard model 
uncertainty is much less in the energy regime probed by 
this electric field sensitivity. Finally, we show the fre- 



quency dependence of eq. ( 11 ) using the standard model 



cross section for neutrinos and the high cosmogenic neu- 
trino flux. In Fig. |9j we show the predicted number 
of events for 100 hours as a function of detected radio 
frequency. For standard model neutrino nucleon cross 
sections and the cosmogenic neutrino flux, the electric 



IV. LUNAR RADIO CHERENKOV SIGNALS 
FROM COSMIC RAY PROTONS 

In the previous section, we saw that for a range of 
energies, depending on detection parameters, the down- 
ward neutrinos dominate the event rate. Contributions 
in the smooth case can dominate for some energies and 
for lower radio frequencies, while the surface roughness 
is important for the higher radio frequencies. Given that 
downward production of hadronic showers by neutrino 
interactions in the lunar regolith can produce observable 
radio Cherenkov signals at Earth, one should also con- 
sider the corresponding signals from hadronic showers 
induced by cosmic rays, the topic of this section. The 
Westebork group has already used the absence of a cos- 
mic ray induced radio Cherenkov signal to put a limit on 
the cosmic ray flux |29| . 

Our starting point is to treat the cosmic rays incident 
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FIG. 7: The number of events for a 1.5 GHz signal with £ m i n = 
10 -8 V/m/MHz with the high cosmogenic flux from Fig. I] 
with E v < 10 14 GeV, as a function of S = a vN /(JuN '• The 
up, down smooth and down rough contributions are as in Fig. 
[5] (b). The dots labeled CTW show the uncertainty bands 
ofRef. M- 
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FIG. 8: The number of events for a 1.5 GHz signal as in 
Fig. [7} with e mill = 10" 11 V/m/MHz with the low and high 
cosmogenic flux from Fig. 
function of S. 



with E v < 10 14 GeV, 



on the lunar regolith with a flux shown in Fig. [T] The 
Pierre Auger Cosmic Ray Observatory analysis favors a 
composition of primarily iron nuclei at ultrahigh energies 
|38j . Iron nuclei incident with the same cosmic ray energy 
as a single proton will produce similar hadronic showers, 
so we use the cosmic ray flux of Ref. [3U] assuming that 
the incident particles are protons carrying all the energy. 
For neutrino induced hadronic showers, we approximate 
the hadronic shower energy to be E s } lr ~ 0.2E„. We will 
make the same approximation for cosmic rays, E s h r ~ 
0.2E CR . 

The essential difference between incident cosmic rays 



FIG. 9: The number of neutrino events as a function of radio 
frequency from cosmogenic neutrinos (high flux) with stan- 



dard model interactions for e n 



10" 



1(T S V/m/MHz. 



and neutrinos is the difference between strong interaction 
and weak interaction cross sections. For E v = 10 12 GeV, 
the neutrino interaction length in the standard model is 
on the order of 10 7 cmwe, as indicated in Fig. 2, how- 
ever, for cosmic rays, the interaction length is Lqr — 50 
cmwe. This new scale makes some of the analytic ap- 
proximations in Ref. [27] inapplicable to cosmic rays. 

The short interaction length of cosmic rays makes at- 
tenuation of the flux important for downward cosmic 
rays, and it completely extinguishes the upward flux. We 
consider here the modifications to the effective aperture 
for the downward flux contribution, including attenua- 
tion in the regolith. Details of our evaluation appear in 
Appendix A. 

Without attenuation, the target volume integral is gov- 
erned by the maximum depth related to the photon at- 
tenuation length L~ t . In Ref. [27], the integral over the 
diameter of the moon r is replaced by an integral over the 
perpendicular distance to the surface h. In the small an- 
gle approximation, the maximum depth is approximately 



hmax — L~f sin# c /Q I 



A 2 



f 2 A 2 



(19) 



where A is the polar angle from the Cherenkov peak. (See 
eq. (13) in Ref. [27].) This factor of L 7 arises because of 
the requirement that the radio Cherenkov signal emerge 
from the regolith. Without flux attenuation, the integral 
over h simply contributes a factor of h mncr . These are 



some of the factors that appear in eq. ( 12 1 for neutrinos. 



With flux attenuation, the integral over h in the eval- 
uation of the effective aperture results in a factor of 
I sin a I Lcr, where a is the angle of the incident cosmic 
ray with respect to the horizontal, as in Ref. \27 . This 
modification of eq. (29) in Ref. [57], and subsequent 
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approximate integration, yields 
v/^l(/oAo) 3 



Pcr(E) 



12 



1+T- 



(20) 



The first term in parenthesis is the smooth contribution, 
and the second term includes the additional contribution 
from surface roughness. 

Eq. ( 20 ) shows that the probability for a cosmic ray to 



produce a signal is independent of the cosmic ray cross 
section, as long as cosmic ray flux attenuation is impor- 
tant on the scale of L 7 . This comes from two compen- 
sating factors. One factor of Lqr comes from the limit 
on the depth of targets from which there a signal, not 
because of radio wave attenuation, but because cosmic 
rays do not penetrate deeper in the regolith. The second 
factor is proportional to ucr ~ I/Lcr from the prob- 
ability that the cosmic ray interacts to produce a radio 
Cherenkov signal. 

In Fig. [lOj we show the effective aperture for incident 
cosmic rays for v = 0.15 and 1.5 GHz and e m in = 10~ 8 



V/m/MHz, and in Fig. 11 we show the same for e m i n = 
10~ u V/m/MHz. We see a similar dominance of the 
smooth contribution for the lower frequency, and rough 
contribution for the higher frequency, in the cosmic ray 
induced radio Cherenkov signals. 

Based on these results for the effective aperture, we 
evaluate the number of events as a function of radio fre- 
quency induced by cosmic ray interactions. They are 
shown in Fig. [12] with the solid lines. For reference, 
the neutrino induced Cherenkov rates are shown with 
dashed lines. In all cases, the cosmic ray induced rates 
are larger than the neutrino induced rates. This is true at 
low and high radio frequencies, where either the smooth 
or rough downward contributions dominate. Our conclu- 
sion is that in the standard model of neutrino interactions 
with nucleons, the cosmic ray induced Cherenkov signals 
will overwhelm the cosmogenic neutrino induced signals, 
at least for the sample fluxes shown here. The lower 
electric field thresholds sample lower neutrino energies 
where there are larger fluxes of neutrinos, but there are 
even larger fluxes of cosmic rays at those energies. The 
higher energy thresholds suffer from low rates due to the 
low flux. 

Theoretical predictions for neutrino induced event 
rates are enhanced if either the neutrino flux is much 
larger than shown, for example, in Fig. [T] or if the 
neutrino cross section is much larger than the standard 
model cross section. It is this second possibility that we 
explore in the next section. 



V. RADIO CHERENKOV SIGNALS AND 
ENHANCED NEUTRINO CROSS SECTIONS 

Enhanced neutrino-nucleon cross sections can arise in 
a variety of extensions of the standard model. Here, we 
consider the case of large extra dimensions. In theoretical 
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FIG. 10: The effective aperture as a function of cosmic ray 
energy for (e min = 10~ 8 V/m/MHz for (a) v = 0.15 GHz 
and (h)u = 1.5 GHz. The solid line shows the total of the 
downward rough (dr, dot-dashed) and downward smooth (ds, 
dotted). The cosmic ray effective aperture is independent of 
the cosmic ray cross section. 



models with large extra dimensions and low scale grav- 
ity there is a possibility of creating a microscopic black 
hole in neutrino-nucleon interactions at very high ener- 
gies [33]. In these models, gravitational interactions are 
modified and the four dimensional Planck scale (Mpi) is 
related to the fundamental Planck scale in 4 + No di- 
mension (Md) by 

Afpi = M% d+2 V Nd 

where Vn d — (2itR) Nd is the volume of the 7V£>-torus 
and R is the size of extra dimensions [3j3 0DJ ■ When R 
is large (of the order of a millimeter) , and the number of 
extra dimensions No is larger than 2, the fundamental 
Planck scale can be of the order of few TeV. 

In high energy collisions, when particles with energies 
above Mo approach each other at the impact parameter 
which is less than the Schwarzschild radius in 4 + No 
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FIG. 11: The effective aperture as a function of neutrino en- 
ergy for (e min = KT 11 V/m/MHz for (a) v = 0.15 GHz and 
(b)u = 1.5 GHz. The solid line shows the total of the down- 
ward rough (dr, dot-dashed) and downward smooth (ds, dot- 
ted). 



dimensions, they can form the 4 + Nd dimensional black 
hole with mass Mbh [HI 14TH45] . The radius rg is given 
by 06] 



1 1 



rs 



fjrM 



D 



M BH f8T(^y 



Mr 



Nr 



(21) 



The neutrino-nucleon cross section for black hole produc- 
tion is given by [32] 

a{vN -> BH) = J2 f {MTS) , dx af H (xs) f t (x, Q 2 ), 

(22) 

where af H is the neutrino-parton cross section given by 



s is the center of mass energy squared, s — 2mjqE v , and 
fi(x,Q 2 ) is the parton distribution function for parton 
i. For semiclassical approximation to be valid, we need 
Mgjf » M D . 

Current limits on the Md and Njy come from collider 
data [35J 27] as well as from the astrophysical observa- 
tions. Strongest limits for Nd < 4 come from super- 
nova cooling and neutron star heating, Mjj > 4 TeV for 
N D = 4, M D > 0.8 TeV for N D = 5 and for lower val- 
ues of Nd, Mo is constrained to be much larger than 
few hundred TeV. Non-observation of the black hole pro- 
duction in cosmic neutrinos provide stringent limit for 
Nd > 5, Mo > 1 TeV [25]. 

We consider the parameter space (M D , N D ), for fixed 
Mg"jP using the cross sections from Ref. [HJ where 
M-gjP = Md , to illustrate the consequences of enhanced 
neutrino nucleon cross sections given the highest cosmo- 
genic neutrino flux shown in Fig. 1. 

The neutrino-nucleon cross section for black hole pro- 
duction exceeds standard model cross section for neu- 
trino energies above 10 6 GeV and is about two orders 
of magnitude larger than the st and ard model cross sec- 
tion at E v - 10 11 GeV. In Fig. [l3|we show the effect of 
the mini-black hole contribution to the neutrino-nucleon 
cross section on neutrino interaction length. The dashed 
line shows the interaction length for No = 7 and Md = 1 
TeV as calculated in Ref [13] . The dotted line shows 
the interaction length for the same (Nd, Md) with the 
cross section taken from Ref. [U] . We note that the neu- 
trino interaction length due to black hole production had 
strong energy dependence and at neutrino energies above 
10 9 GeV, it is order of magnitude smaller than in the case 
of the standard model neutrino interactions. We use the 
neutrino cross sections due to black hole production from 
Ref. [13] as sample cross sections in what follows below. 

The fact that the neutrino interaction length decreases 
to less than 10 4 cmwe above E v — 10 12 GeV means that 
neither the cosmic ray approximation of strong attenua- 
tion nor the standard model neutrino evaluation of the 
effective downward aperture applies. The details of the 
evaluation of the effective aperture with the mini-black 
hole enhanced neutrino cross section contribution to the 
probability function Pbhu(E) is shown in the Appendix. 

To illustrate dependence of Pbhu(E) on the neutrino 
interaction length, we evaluate Pbhv {E) a s a function of 
the ratio of the L 1 / L v 



We show in Fig. 14 with the solid 



a{uj ->■ BH) = tttUMbh = Vs)6(Vs - M%g) , (23) 



line the numerical evaluation of Pbhu(E) from the for- 
mulas in GMJ with the modification of neutrino attenua- 
tion even over the distance scale of L 7 . The dashed hori- 
zontal lines show the what could be called the "cosmic ray 
limit, " when the non-standard model interactions make 
the neutrino interaction length small compared to L 7 . 
The dot-dashed lines show the "standard model neutrino 
limit" when L v ^> L~ f (the analytic expression of GMJ). 
In order from top to bottom, the curves show (e m i n) v) for 
(KT 11 , 0.15), (10~ 8 , 0.15), (HT 11 , 1.5) and (10~ 8 , 1.5) in 
units of (V/m/MHz, GHz). We note that the probability 
function for enhanced neutrino cross sections approaches 
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the cosmic ray approximation for L 7 / ' L v > 0.1 (1.0) when 
£min = 10" 8 V/m/MHz (e min = IO" 11 V/m/MHz) and 
v = 1.5 GHz (i/ = 0.15 GHz). For lower values of L~/L v , 
the probability follows dependence on L 7 /L„ similar to 
the the downward neutrino approximation of Eq. (12). 
As noted in Ref. [27], the approximate analytic solu- 
tions reproduced in our Eq. ( 12 ) are good to better than 



25% compared to a numerical solution of the multidimen- 
sional integrals in their analytic expression. In evaluating 
neutrino induced Cherenkov signals with mini-black hole 
enhanced cross sections, we will use the numerical result 
(the solid line in Fig. 14) to ensure that we are prop- 
erly accounting for the transition from weak to strong 
neutrino cross sections. 
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FIG. 12: The number of events as a function of radio fre- 
quency from the high flux of cosmogenic neutrinos (dashed) 
and from cosmic rays (solid), from top to bottom with e m i n = 
10~ n — 10 -8 V/m/MHz for standard model neutrino inter- 
actions. 



FIG. 13: The neutrino interaction length, with contributions 
from non-standard model mini-black hole production from 
Connolly, Thorne and Waters [H| and Feng gT] for N D = 7, 
M D = 1 TeV. 




In Fig. [15] we show number of cosmogenic neutrino 
events in 100 hrs for two frequency choices, v = 0.15 
GHz (dashed) and v = 1.5 GHz (solid), as a function of 
the minimum electric field £ m i n . The parameters {Np, 
M D ) are (1, 1 TeV), (7, 2 TeV) and (7, 1 TeV) from low- 
est to highest. The highest pair of solid and dashed lines 
for low e m i n is for the cosmic ray induced events. We find 
that for v = 1.5 GHz (solid lines), the neutrino induced 
BH contributions are larger than the cosmic ray contri- 
butions for e min > 2 x 10~ 9 V/m/MHz (e min > 6 x 10~ 9 
V/m/MHz) for N D = 7 and and M D = 1 TeV (M D = 2 
TeV) black hole parameters. This comes from the high- 
est energies considered, where the neutrino nucleon cross 
section is much bigger than the standard model, so the 
standard model neutrino interactions are negligible. 

For v = 0.15 GHz, shown with the solid lines in Fig. 
[15] the neutrino induced BH contributions are larger than 
the cosmic ray induced Cherenkov signal for e m j n > 10~ lf} 
V/m/MHz, lower than the £ m i n for v = 1.5 GHz required 
for a signal from neutrinos with Np, = 7 and Mp. = 1 
TeV black hole parameters. 



FIG. 14: The probability function which enters into the ef- 
fective aperture, for neutrinos in which there is an enhance- 
ment of the neutrino cross section. The solid line shows the 
numerical evaluation, assuming E = 10 14 GeV, but varying 
the neutrino cross section relative to a fixed L 7 . The dashed 
horizontal line is the cosmic ray approximation of Eq. ( |20[ ) 
and the dot-dashed line is the downward neutrino approxi- 
mation of Eq. ( |12[ ). In order from top to bottom, the curves 
show (emin,!') for (10 -11 , 0.15), (IO -8 , 0.15), (10 -11 , 1.5) and 
(IO -8 , 1.5) in units of (V/m/MHz, GHz). 



From Fig. 15 



we note that different BH parameters 
give similar dependence of the events on minimum elec- 
tric field, with only difference being in the overall number 
of events. With 100 hrs, the BH enhanced neutrino sig- 
nals dominate the cosmic ray background in the range 
of 40 events at the crossover for v = 0.15 GHz, while 
for v = 1.5 GHz, the crossover occurs at a fraction of 
an event level for the largest BH cross section considered 
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16 we show the frequency dependence of the 



here. 
In Fig. 

event rate for the lower and higher minimum electric field 
for three values of (Np, Mr,)- There is roughly an or- 
der of magnitude change in the neutrino induced event 
rates as a function of radio frequency for v = 0.15 — 1.5 
GHz. Again we see in the figure that while event rates 
are higher for lower minimum electric fields, the cos- 
mic ray background is much larger than the signal and 
thus dramatically lowering minimum electric fields does 
not favor the observation of cosmogenic neutrino induced 
Cherenkov signals. 

Our conclusion is that for cosmogenic neutrino fluxes 
and neutrino cross sections enhanced by BH production 
with Nn — 7 and Mr, = 1 TcV, the neutrino signal equals 
the cosmic ray background for £ m ; n ~ 1CT 10 V/m/MHz 
when v = 0.15 GHz, and the number of events is about 
40. Thus, improvements in detector sensistivies would 
be required to observe the neutrino black hole enhanced 
signal. For a current detector capability of e min ~ 10~ 8 
V/m/MHz, the enhanced neutrino cross section domi- 
nates the radio Cherenkov signal from the Moon relative 
to the cosmic ray background, however, the number of 
events in 100 hours is too small, of the order of 4 x 10~ 4 . 

In addition to the cosmogenic neutrino flux that we 
have considered, there is a possibility that astrophysical 
sources produce larger neutrino fluxes via Fermi shock ac- 
celeration of the charged particles such as protons which 
collide with other protons or photons in a disk or a jet 
producing mesons which decay into neutrinos. Neutrino 
fluxes due to shock acceleration follow a power-law, i.e., 
~ E~ 2 . We consider a neutrino flux which is cur- 
rently below the Anita limit [TB], i.e. , 



E 2 § v = 5 x 10~ 7 GeV/cm 2 /s/sr 



(24) 



In Fig. 17 we show the event rates for the neutrino flux 



in eq. ( |24[ ) and for neutrino enhanced cross sections due 
to black hole production obtained with different black 
hole parameters (Nr>, Mr,) from Connolly et al. [TJ], as 
a function of the minimum electric field. We find that 
even for e m j n ~ 10~ 8 V/m/MHz, which is a current de- 
tector capability, for any frequency above 150 MHz the 
signal is above the cosmic background by several orders 
of magnitude. Total number of events varies between 1 
(for N D = 1 and M D = 1 TeV) and 20 (for N D = 7 and 
M D = 1 TeV) for v = 1.5 GHz and between 40 and 200 
for v — 150 MHz, depending on BH parameters. Low- 
ering e m ; n results in only a slight increase in the event 
rates. For £ m ; n > 4 x 10~ 10 the neutrino signal is above 
the cosmic ray background for any choice of BH param- 
eters that we considered. 



VI. DISCUSSION 

We have used the formalism of GMJ in Ref. [37] 
to evaluate the effective aperture for lunar observations 



from Earth of radio Cherenkov signals produced by 
neutrino interactions. The GMJ analytic effective aper- 
ture agrees qualitatively with the Monte Carlo results of 
James and Protheroe [15], as discussed in Ref. [57]. In 
the analytic approach, we have made a new evaluation of 
the lunar effective aperture for cosmic rays. This relics 
on the result of Ref. [35] that cosmic ray interactions in 
the lunar regolith produce radio Cherenkov signals that 
are indistinguishable from neutrino induced signals. 

As Fig. 12 shows, the cosmic ray induced event rate 
always dominates the standard model neutrino event rate 
for the cosmogenic models of Ref. [9 for e = 10~ 8 — 10~ n 
V/m/MHz. Current capabilities for the electric field 
threshold are on the order of e m i n ~ 10~ 8 V/m/MHz. 
Decreasing £ m ; n will not improve the capability of this 
lunar technique to detect the cosmogenic neutrino flux if 
the standard model interactions, e.g., with cross sections 
on the order of those reported in Refs. [TO] and [T5], are 
the only interactions responsible for the neutrino cross 
section. While lowering the minimum electric field in- 
creases the neutrino induced signal, it also increases the 
cosmic ray induced signal. The effective aperture for cos- 
mic rays is independent of the cosmic ray cross section. 
Thus, the current discussion about the cosmic ray com- 
position at the highest energies 38J will not impact our 
results as long as the cosmic rays interact strongly. 

The radio Cherenkov technique can be exploited for 
lunar observations if the neutrino nucleon cross section is 
increased. We showed that for the choice of Mr, = 1 TeV 
and Nr> = 7 for BH production, an observable neutrino 
induced radio Cherenkov signal from the Moon is induced 
and it is larger than the cosmic ray induced background. 
For v = 150 MHz, this would require an instrumental 
improvement to reach e m i n ~ 10~ 10 V/m/MHz for the 
high cosmogenic flux of Fig. 1. We note that eventually 
enhancements of the neutrino cross section will push the 
rate evaluation into the "cosmic ray regime," where the 
rate is independent of the cross section. 

In Fig. [IT] we showed that if both the neutrino cross 
section and the neutrino flux are enhanced, lunar radio 
Cherenkov techniques can make observations rather t han 
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set limits, even for £ m j n ~ 10 -8 V/m/MHz. Fig. 
shows the observational capability as a function of e 
for two frequencies at the current limit on a neutrino flux 
of E 2 $ v = 5 x 10~ 7 GeV/cm 2 /s/sr for a specific mini- 
blackhole model of neutrino cross section. In this model, 
the cross section enhancement relative to the standard 



model is S = o w n / c 



SM 
uN 



500 for E v = 10 12 GcV. 



Taking a more generic approach, looking at the de- 
tectability as a function of S, we can see the capabilities 
of the lunar technique based on the GMJ formalism. In 
order to study the sensitivity of radio Cherenkov detec- 
tion to the neutrino flux, we take neutrino flux of the 
form, 



E 2 <$> V (E V ) =Ax IQ- 8 GeV/cm 2 /s/sr. 



(25) 



We investigate the value of A such that when A = A 



100, 



13 




10" 
10 4 

io 2 

10° 

io 2 
io 4 
10"* 
10 s 



! 

f^iC^i : 

t~ — H^^?^ — ■ 
i ■ ■ 


v = 1 50 MHz : 

v= 1.5 GHz 1 

-^iiiii" - - = - - : : : ~ : : - ~ - - — 

"-^^-i zzzzzizzzir — ^^zzzz~ 




K — — zzzm~~ 














r N = l.M = 1 

! D D 




j- N D =7,M D =2 




r — N D =7,M D =1 




r cr 





10 



io 1 



10 



10 



FIG. 15: Events in 100 hrs for v = 150 MHz (dashed) and 
1.5 GHz (solid) as a function of e m i n (in V/m/MHz) using 
enhanced neutrino cross sections from Connolly et al. [14] , 
for three choices of black hole parameters, Nd and Mo- The 
high cosmogenic flux from Fig. 1 is used. The parameters 
{N D , M D ) are (1, 1 TeV), (7, 2 TeV) and (7, 1 TeV) from 
lowest to highest event rates respectively. The top solid and 
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FIG. 16: Events in 100 hrs as a function of radio frequency 
for £ m i n = 10 -11 V/m/MHz (upper curves) and e m - ln = 10~ 8 
V/m/MHz (lower curves) for three choices of mini-blackhole 
parameters with the high cosmogenic neutrino flux. The cos- 
mic ray induced event rate is shown with the dotted lines. 



the theoretical prediction is one neutrino event for 100 
hours of observation. We consider not only the total 
event rates but also the ratio of the cosmic ray induced 
signal over the neutrino induced signal, 



R 



FCR 



FIG. 17: Events in 100 hrs as a function of e m i n for the 
enhanced neutrino cross section due to black hole produc- 
tion, for different values of (No, Mo) as in Fig. 15 and 
£ 2 $„ = 5x 10" 7 GeV/cm 2 /s/sr. The lowest solid and dashed 
curves at e m i n = 10~ 8 V/m/MHz are the cosmic ray rates. 



In Figs. 18 and 19 we show the two quantities, R and 
the minimum A = ^4ioo req uire d for one neutrino event in 
100 hrs, for the flux in eq. ((251) for e min = 10~ 8 and 10~ 9 



V/m/MHz, respectively, as a function of S = Pun/p^n 
The solid lines in these figures show the ratio R using 
the expression in eq. (A8) which produces the solid lines. 
This is labeled the "numerical" result. The dashed lines 
show the result using the "neutrino approximation" in 
which attenuation of the downward neutrino flux is not 
included. The dot-dashed lines show the value of Aiqq 
required to produce one event in 100 hours. 

For e = 10 -8 V/m/MHz, for both frequencies shown 
(1.5 GHz and 150 MHz), R < 1. The deviation be- 
tween the solid and dashed lines shows the expected re- 
sult that eventually, the neutrino rate does not increase 
with cross section because of attenuation of the down- 
ward flux in the Moon. The dot-dashed lines showing 
the minimum A for a neutrino flux with an analytic form 



(26) 



following eq. (25) required for one neutrino event in 100 
hrs also show the saturation effect due to attenuation. 
Even with large neutrino cross sections (S 3> 1), given 
100 hrs, the minimum observable A for v = 1.5 GHz and 
e = 10 -8 V/m/MHz is A 100 ~ 1. For values of S closer 
to 1, the minimum observable A is several hundred for 
100 hrs of observation. This range of parameter space is 
already excluded by ANITA |17j . however, not in this en- 
ergy range. The lower frequency of v = 150 MHz probes 
a lower value of A, from A = 1 down to A — 0.2. 

In Fig. [19] we show the same ratio R and A wo for 
e = 10 -9 V/m/MHz. In this case the values of Aioo 
needed to reach one neutrino event are much lower than 
in the case when e = 10 -8 V/m/MHz. This is due to 
the fact that lowering s we are probing neutrino fluxes 
at a lower energy where the flux is higher, as can be seen 
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tween 10MHz and 80 MHz, is expected to lower their 
energy threshold down to 10 11 GeV and with 30 days 
of data taking will be probing the neutrino flux down 
to E 2 $ u ~3x 10~ 9 GeV/cm 2 /s/sr assuming standard 
model neutrino interactions |25j . Cosmic ray background 
for LOFAR is still sufficiently small relative to the neu- 
trino signal. The advantage of LOFAR over RESUN, 
for example, is that it operates at low frequency and 
has longer observation time. However, it is important to 
note that lowering energy threshold down to 10 10 GeV 
(i.e. lowering e min below 10~ 10 V/m/MHz), would result 
in larger neutrino rates but the cosmic ray background 
would become significant. For neutrino fluxes lower than 
£ 2 $„ - 3 x 10~ 9 GeV/cm 2 /s/sr, LOFAR provides an 
excellent probe of physics beyond the Standard Model. 



FIG. 18: For E 2 $ u = A X IO -8 GeV/cm 2 /s/sr, the ratio 
R — Vcr/^v for A = 1 is shown with the solid lines where the 
full expression for the neutrino rate is used, while the dashed 
lines are evaluated using the "neutrino approximation," with 
no attenuation for down-going neutrinos, as a function of S = 
(TuN/crSff for e = 10 -8 Vm -1 MHz -1 . The dot-dashed lines 
show the minimum A to produce one neutrino event in 100 
hrs, as a function of S. 



from Fig. 1. However, for v — 1.5 GHz, the ratio R is 
also much larger, and for S < 10 2 , the CR background 
is larger than the neutrino signal. For lower frequencies, 
when v = 150 MHz, for any value of S, including 5 = 1 
which corresponds to the standard model neutrino cross 
section, the CR background is much smaller than the 
neutrino signal. In general, lower frequencies give larger 
event rates, and operating at lower frequency has this 
advantage in detecting neutrino fluxes. A theoretical ad- 
vantage at the lower frequencies is that the event rate 
relies less on the surface roughness than for the higher 
frequencies. However, the observational challenge of at- 
mospheric dispersion of the signal at v = 150 MHz is 
much more significant than at v = 1.5 GHz. 

There are several radio telescopes that have recently 
looked for very high energy neutrinos using the Moon 
as a target. The RESUN Project used the extended 
very large array (EVLA) with v ~ 1.5 GHz and e m i n = 
10~ 8 V/m/MHz [3T]- The Westebork Synthesis Ra- 
dio Telescope (WSRT) operates in the frequency range 
of 115 — 180 MHz and has reported limit on neutrino 
flux of E 2 <$> v ~ 10~ 6 GcV/cm 2 /s/sr [22 . The Lunaska 
experiment at the Australia Telescope Compact Array 
(ATCA) [23] is operating in the region v = 1.2 - 1.8 GHz 
and is expected to have sensitivities to lower neutrino 
fluxes than RESUN. At ultra-high energies [E v > 10 14 
GeV), NuMoon [53] which uses the Westebork Synthe- 
sis radio Telescope in Netherlands, one of the most sen- 
sitive low-frequency experiment (v — 113 — 175 MHz), 
will be able to place more stringent limit on the neu- 
trino flux. LOFAR (Low Frequency Array), which cov- 
ers frequencies between 120 MHz and 240 MHz and be- 
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FIG. 19: The ratio R and the minimum Axoo, as in Fig. [18] 
but for s = 10 -9 Vm _1 MHz~ 1 . 
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Appendix A: Effective aperture for downward 
particle fluxes 

Following Gayley, Jaeger and Mutel in Ref. [37], we 
evaluate the probability function P{E) for downward in- 
cident particles with 

1 L f° r°° 
P{E) = --2 / dacosa / dA sm(6 c + A) 

77 Lv J— oo i-oo 

" dze- T "H R H D i , (Al) 
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(eq. (12) of Ref. [17]). For the sake of discussion, we will 
call the incident particles neutrinos. Here, a is the angle 
the incident neutrino makes with respect to the horizon- 
tal, with a < indicating "downward" neutrinos. The 
quantity A characterizes the polar angle of the interior 
ray solid angle relative to the cherenkov angle 6 C of the 
radio signal, and <j> is its azimuthal angle. The factor £ 
accounts for the fact that the rado signal encounters the 
lunar surface from the inside. The quantity z is the depth 
h in the lunar regolith, normalized by L~ r In this equa- 
tion, % r and T-L d are there to satisfy the conditions that 
the radio ray refracts (and does not totally internally 
reflect inside the moon), and that the rays are bright 
enough (detectable). Ref. [27] has Hr written in terms 
of two more integrals, accounting for surface roughness 
so that the surface tilt polar angle is a = <jqw and the 
azimuthal tilt direction is 4>' relative to a smooth surface, 
with 



2 



tt/2 



dw 



(A2) 



The Heaviside step function r hi WiV hi i enforces the require- 
ment that the ray emerge from inside the Moon. Finally, 
the factor e~ T " accounts for the attenuation of the inci- 
dent neutrino flux with 

Tjj = z • L 1 jL v ■ 1/| sina| . (A3) 

The integral over z is straightforward, with 



dze -r„ = I sina|L^ / 1 c _ 2 , ( .„ hl „ L , , 



L, 



The quantity z max is 



sin0 c / o 2 1 



A 2 



f 2 A 2 



(A4) 
(A5) 



The result in eq. (A4) is the origin of the two limits, 



the "cosmic ray limit" and the "standard model neutrino 



limit." In the cosmic ray limit, 
Ismail,,, 

dze " ~ ■ 

/o 



/' 

Jo 



CR limit 



(A6) 



since the argument of the exponential is a large negative 
number with L 7 /L„ large. The neutrino limit involves 
the expansion of the exponential in a power series, keep- 
ing the first non-zero term, so 



dze 



neutrino limit 



(A7) 



For enhanced neutrino cross sections, where neither 
the cosmic ray nor neutrino limits are applicable, we nu- 



merically integrate eq. ( Al ) in the small angle limit, with 
the Heaviside functions enforced with downward incident 
angles to give 
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where a m i n = — A — wctqCos^' — — A — w/wq. Ana- 
lytically, in the downward neutrino limit with sin#o = 
\l"n%. — l/n r , this yields 



n r — 1 L 7 3 A /it a 16 



P ^=^^-^^oAo(/oA + ^o 
In the cosmic ray limit, the integrals give 



P(E) = P C r(E) 



y^U/oAo) 3 / 3 al 
12 I + 4/ 2 A 2 
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